AD**A072  548 


UNCLASSIFIED 


CALIFORNIA  STATE  UNIV  LONG  BEACH  DEPT  OF  MECHANICAL  ->ETC  F/6  20/4 
PREDICTION  OF  UNSTEADY  TURBULENT  BOUNDARY  LAYERS  WITH  FLOW  REVE— ETC(U) 
1979  T CEBECI*  L W CARR>  P BRADSHAW  NOOO  19-78*^0466 

NL 


1 


PREDICTION  OF  ^NSTEADY  TURBULENT  BOUNDARY^YERS  / 
/ WITH  , FLOW  REVERSAL/-)  / 


Tuncer  /^ebeci' 


Ca  90840 


ffett  Field,  Ca  94035 


Presented  at  the  Second  Symposium  on  Turbulent  Shear  Flows 
July  2-4,  1979,  Imperial  College,  London,  England 


T^i  documast  has 
for  |»ukBe  < 

awrftfUai  It  vaMi 


08  06 


I 


I 


PRtDICTlON  Of  IINSTIADY  TUKHULtNT  BOONUARV  LAYLRS  WITH  FLOW  REVERSAL* 


Tuncor  Ccbeci  / 

Professor  of  Mechanical  Imjineerinq , California  State  University  at  Loikj  Beach,  Ca.  90HA0 

Lawrence  W.  Larr 

Reseirch  Scientist,  U.S.  Amiy,  AeronH'rhan ics  lab.,  Moffett  Meld,  Ca.  9403S 

Peter  Bradshaw 

Professor  of  Experiniental  Aerodynamics,  Imperial  College,  London,  England 


INTRODUCTION 

%)The  prediction  of  unsteady  turbulent  boundary 
layers  with  flow  reversal  is  of  importance  in  a num- 
ber of  aerodynamic  problems,  notably  in  dynamic  stall, 
buffeting  and  gust  studies.  However,  some  of  the  more 
popular  turbulence  models  implicitly  assume  that  the 
wall  shear  is  positive  and  their  extension  to  un- 
steady flows  with  flow  reversal  is  not  easy.  It 
requires  a modification  of  the  functional  form  of  the 
law  of  the  wall  and  the  manner  in  wiiich  ttie  wall 
shear  is  determined.  Two  near-wall  assumptions  are 
considered  here.  In  the  first,  the  near  wall  grid 
point  is  located  in  the  logarithmic  region  and  the 
law  of  the  wall  is  used  to  link  the  flow  properties 
at  this  grid  point  to  the  wall.  In  the  second,  a , 

Van  Driest  formulation  due  to  Cebei  i and  Smith  (it" 
is  used  and  implies  that  the  grid  point  closest  to 
the  wall  will  occur  in  the  viscous  sublayer. 

A further  aspect  of  these  flows  o(  current^ nter- 
est  is  the  possibility  of  a singularity  occurring  in 
the  reversed-flow  region.  Examples  of  this  phenom- 
enon have  also  been  reported  in  laminar  flows  but, 
in  earlier  studies  two  of  us  (Cebeci  (2^,^)  and 
Bradshaw  (4))  have  shown  that  the  occurrence  is  not 
a feature  of  the  governing  equations  hut  is  due  to 
the  limitations  of  the  numerical  procedure  used.  We 
shall  demonstrate  that,  for  the  examples  we  study, 
there  is  no  indicatation  of  such  a singularity  in 
turbulent  flow  either  but  there  is  a clear  indication 
of  rapid  thickening  of  the  boundary  layer. 

In  addition  to  the  examination  of  wall  functions, 
we  have  also  considered  two  turbulence  models  for 
unsteady  flows  without  flow  reversal.  The  algebraic 
eddy-viscosity  formulation  of  Cebeci  and  Smith  (CS) 
is  compared  with  the  transport  model  of  Bradshaw, 
Ferriss  and  Atwell  (5)  (BF).  Calculations  were  per- 
formed to  determine  whether  the  representation  of 
unsteady  flows  with  strong  pressure  gradients  requires 
that  account  be  taken  of  transport  of  turbulence 
quantities.  As  will  be  shown,  the  predictions  with 
both  models  are  nearly  identical  for  steady  flows  but 
indicate  significant  differences  with  the  increasing 
pressure  gradient  for  unsteady  flows, 

In  the  following  sections  we  describe  the  turbu- 
lence modelling,  the  governing  equations,  the  numer- 
ical procedure  and  the  results,  respectively. 

TURBULENCE  MODELLING 

The  Reynolds  stresses  are  governed  by  the  well- 
known  transport  equations  whose  left-hand  sides  are 


"substantial"  derivatives  of  the  Reynolds  stresses  - 
that  is,  the  rate  of  change  of  Reynolds  stress  with 
respect  to  time  as  seen  by  an  observer  following  the 
mean  motion  of  the  fluid  substance.  We  write,  for 
example 
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where  in  the  case  of  an  unsteady  flow  u,v  and  w 
are  defined  as  ensemble  averages  over  many  repeats  of 
the  same  experiment  and  are  functions  of  time  measured 
from  the  instant  at  which  the  unsteady  flow  is  started. 
Here  u'  and  v'  arc  fluctuations  about  the  ensemble- 
averages  u and  v,  and  are  quasi-random  functions 
of  time  (and  spatial  position)  as  usual  , while  u'v' 
is  Itself  an  ensemble  average  and  thus  a relatively 
smooth  function  of  time.  The  riqht-hand  sides  of  the 
Reynolds  stress  transport  equations  are  exactly  the 
sanK'  as  in  steady  flow  if  the  overbars  are  interpreted 
as  ensemble  averages.  Therefore,  the  relative  sue 
of  the  time-derivative  and  space-derivative  contribu- 
tions to  the  left-hand  side  of  the  Reynolds  stress 
transport  equations  is  immaterial  to  a turbulence 
model  based  on  term-by-teTiii  empirical  approximation  to 
the  right-hand  side  o*  the  Reynolds  stress  transport 
equations.  All  that  matters  is  the  size  of  Du'v'/Ot 
(hereafter  referred  to  simply  as  "the  rate  of  change 
of  shear  stress")  as  a whole.  Therefore,  the  ability 
of  a turbulence  model  to  predict  unsteady  flows  can  be 
gauged  by  its  ability  to  predict  steady  flows;  com- 
monly, turbulence  models  become  less  satisfactory  when 
the  rate  of  change  of  shear  stress  is  large,  because 
the  empirical  approximations  on  the  riqht-hand  sides 
are  implicitly  based  on  assumptions  that  the  turbu- 
lence structure  parameters  change  only  slowly.  Some 
calculation  methods  are  implicitly  based  on  the 
assumption  that  the  rate  of  change  of  shear  stress  is 
itself  small;  examples  include  the  eddy-viscosity 
assumptions  used  in  some  of  the  calculations  presented 
below.  Again,  however,  the  potential  of  the  method 
for  predicting  unsteady  flows  can  be  gauged  from  its 
ability  to  predict  steady  flows  with  rapid  changes  in 
the  streamwise  direction. 

It  is  not  possible  to  grade  turbulence  models  for 
unsteady  flow  simply  by  the  largest  rate  of  change  of 
shear  stress  - compared,  say,  to  the  major  term 
v*^  3u/ay  on  the  riqht-hand  side  of  the  equation  - 
that  they  will  reliably  predict.  This  is  because  the 
frcqueocy,  as  well  as  the  magnitude  of  the  time- 
dependent  terms  is  important.  Again,  however,  a 
unified  analysis  for  steady  and  unsteady  flow  is 
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possible;  one  need  merely  consider  the  ratio  of  the 
frequency  seen  by  an  observer  following  the  mean 
motion  of  the  fluid  to  the  typical  eddy  frequency. 
First,  consider  steady  flow  with  — say  — streamwise 
periodicity  of  wavelength  A and  typical  mean  speed 
U;  the  typical  turbulence  frequency  could  be  written 
as  (-U  V )Vl.,  where  L is  the  length  scale,  and 
is  of  the  same  order  as  3u/3y  (equal  to  it  if  L 
happens  to  be  equal  to  the  mixing  length  t),  while 
3u/3y  in  the  outer  region  of  a boundary  layer  is  of 
order  U/4.  The  mean-flow  frequency  seen  by  the  mov- 
ing observer  is  u/A  and  the  ratio  of  this  to  the 
typi.al  turbulence  frequency  is  simply  a/a.  There- 
fore, if  a turbulence  model  will  adequately  predict 
a steady  streamwise-periodic  flow  with  a given  value 
of  4/A,  r say,  the  highest  frequency  of  time- 
dependent,  spatial ly-independent  flow  for  which  the 
model  can  be  tested  is  r times  the  turbulence  fre- 
quency (i.e.  r U/s).  One  would  expect  r to  be  of 
order  0.1  at  most.  In  a time-dependent  flow  with 
rapid  spatial  changes  in  the  streamwise  direction, 
the  combined  moving-axis  frequency  must  be  evaluated 
to  apply  the  above  criterion;  note  that  in  the  import- 
ant practical  case  of  travelling-wave  flows,  the 
moving-axis  frequency  is  smal 1 . 

In  sunmary,  if  both  the  rate  of  change  of  shear 
stress  and  the  moving-axis  frequency  are  no  longer 
than  the  values  the  turbulence  model  can  handle  in 
steady  flow,  the  model  can  be  used  with  confidence 
and  without  change.  No  useful  advice  can  currently 
be  offered  about  turbulence  models  for  unsteady  flows 
with  rapid  changes.  However,  if  the  rapid  changes 
are  driven  by  strong  pressure  gradients,  the  Reynolds- 
stress  gradients  in  the  outer  part  of  the  boundary 
layer  need  not  be  calculated  very  accurately,  while 
the  inner  layer  is  more  tolerant  of  rapid  changes 
because  the  typical  turbulence  frequency  is  high.  It 
is  unlikely  that  the  “rapid  distortion  approximation" 
will  be  usable  except  in  flows  with  extremely  large 
time  rates  of  change. 

GOVERNING  EQUATIONS 


an  algebraic  eddy-viscosity  formulation  developed 
(for  steady  flows)  by  Cebeci  and  Smith  and  the  other 
is  a transport-equation  model  developed  by  Bradshaw, 
Ferriss  and  Atwell.  In  the  CS  model,  we  write  Eq.  (A) 
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with  two  separate  formulas  for  e In  the  so-called 
inner  region  of  the  boundary  layer  (tm),-  is  defined 
by  the  following  formula: 
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In  the  outer  region  is  defined  by  the  following 
formula  m ^ f 
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The  inner  and  outer  regions  are  established  by  the 
continuity  of  the  eddy-viscosity  formulas. 

In  the  BF  model,  which  is  solved  only  outside 
the  viscous  sublayer,  we  assume  t * -u'v'  and  write 
a single  first-order  partial-differential  equation  for 
it;  the  equation  was  originally  developed  from  the 
turbulent  energy  equation  but  can  be  equally  well 
regarded  as  an  empirical  closure  of  the  exact  shear- 
stress  transport  equation.  It  is 


The  continuity  and  momentum  equations  can  be 
written  for  two-dimensional  unsteady  incompressible 
laminar  or  turbulent  thin  shear  layers  as: 
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(2)  Here  a^  is  a dimensionless  quantity,  V is  a 
velocity  and  L is  the  dissipation  lengtFi  parameter 
specified  algebraically,  as  L/4  • f(n)  with  n • y/4 
and  f(n)  gives  as  an  analytic  fit  to  an  empirical 

(3)  curve,  by 


Mere 
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and  we  recall  that  u*  and  v'  denote  fluctuations 
about  the  ensemble-average  velocity;  u'  and  v‘ 
are  zero  in  unsteady  laminar  flow,  and  v 3u/3y  is 
negligible  outside  the  viscous  sublayer  in  a turbu- 
lent flow.  These  equations  are  subject  to  the  usual 
boundary  conditions,  which  in  the  case  of  boundary 
layers  are 

y-0  u-v-0;  y-.4  u - U^(x,t)  (5) 

The  presence  of  Reynolds  stress  term,  -uV  intro- 
duces an  additional  unknown  to  the  system  given  by 
Eqs.  (2)  to  (5).  In  this  paper  we  present  calcula- 
tions using  two  different  turbulence  models.  One  is 


f(n) 


0.4n  n < 0.18 

0.095  - 0.05b(2n  - 1)^  0.18  < n < 1.1 
0.016  exp(-10(n  - 1.1)1  n i 1.1  (11) 


In  a more  advanced  version  of  this  turbulence  model 
(6)  L itself  is  determined  from  a transport  equation; 
we  are  now  programming  this  version  of  the  model  for 
unsteady  flow. 

The  turbulent  transport  velocity  V.^,  nominally 
(p'u*  ♦ u'v'^'l/u’v'  is  proportional  to  a velocity 
scale  of  the  large  eddies  and  is  chosen  to  be 
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where  g(n)  is 


14.24 


i 


I 


33.3n^{0.184  ♦ 0.832n) 


33.3n^(0.368  ♦ 2.496n^) 


18. 7i,  ♦ 14.60 


n < 0.5 
0.5  ^ n < 1.0 

n > 1 .0 


(13) 


In  the  BF  model  equations,  the  Inner  boundary  condi- 
tions for  (2).  (3)  and  (10)  are  applied  outside  the 
viscous  sublayer,  usually  at  yi  * SOv/u^.  In  the 
steady-flow  study  reported  in  (7),  these  boundary 


are: 
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is  evaluated  from  the  continuity 

equation 

(2),  and  a*  is  evaluated  from  (2)  and  (3),  on  the 


assumption  that  the  velocity  u is  given  by 


(17) 
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Integrating  we  can  write 
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We  can  also  write  (19a)  as 
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or  as 
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where  F « u/u,  at  y • yi  and  a*  comes  from  the 
last  term  in  (19b)  and  is  the  same  as  in  steady  flow. 
Equation  (19c)  now  replaces  (16). 


SOLUTION  PROCEDURE 


for  0 * y < y^;  Eq.  (14)  is,  of  course,  a special 
case  of  (17).  The  evaluation  of  a*  is  discussed  in 
Ref.  (7);  the  last  term  in  (18)  can  be  as  large  as 
half  tFe  second  (pressure-gradient)  term.  In  unsteady 
flow  without  flow  reversal,  we  use  the  same  inner 
"boundary"  conditions  at  yi  * 60v/u,,  but  because  of 
the  presence  of  the  time-dependent  term  in  (3),  a* 
becomes  more  complicated.  If  we  again  assume  that 
(17)  holds  - remember  that  the  turbulence  structure 
of  the  inner  layer  is  unlikely  to  be  affected  unless 
the  external-stream  frequency  is  very  high  - as  (2) 
and  (3)  give 


We  used  Keller's  two-point  finite-difference 
niethod  (called  the  Box  method)  to  solve  the  system  of 
equations  described  in  the  previous  section.  The 
application  of  this  method  to  unsteady  flows  with  no 
flow  reversal  using  the  CS  model  has  been  described  in 
Ref.  (8).  Its  application  to  steady  two-dimensional 
flows  using  the  BF  model  is  described  in  Ref.  (2). 

More  we  present  a brief  description  of  its  extension 
to  unsteady  two-dinx'iisional  turbulent  flows  with  tlow 
reversal.  Due  to  the  space  limitations,  we  shall 
also  limit  the  description  to  the  equations  that  use 
the  CS  model.  In  a subsequent  report,  we  hope  to 
present  a detailed  description  of  the  solution  pro- 
cedure for  both  model  equations. 

As  in  our  previous  studies,  we  transform  the 
equations  with 


t • tUg/L,  n * (Up/v7)*‘y  (20a) 


x • x/L, 

and  a dimensionless  stream  function  f(x,n,t),  where 
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Here  u,,  is  a reference  velocity,  L a reference 
length,  and  i|i  is  the  usual  definition  of  stream  func- 
tion which  satisfies  the  continuity  equation  (2).  With 
the  relations  defined  by  (20)  and  with  the  definition 
of  eddy  viscosity,  equations  (2)  to  (4)  and  the  bound- 
ary conditions  can  be  written  as 


1 


(bf“)  4 ^ ff"  + m3 


m,  • 7(f'  ^ - f 

■ 2I4 

f) 

(21) 

\ 9x 

9x 

3t  ' 

0;  n -►  n„,  f'  ’’ 

■ It 

(22) 

• 0,  f • f 

Primes  denote  differentiation  with  respect  to  n 


and 


f • u/u^. 


^“e  ^‘'e^ 

m,  • xlu,  — ♦ —^1 

^ V * 97  »T  ' 


(23) 


b • 1 4 r , 
m 


C*- 

m V 


For  simplicity,  we  shall  now  drop  the  bars  on  x and  t. 
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standard  Box  Method 

To  solve  (21)  and  (22)  by  the  standard  Box 
method,  we  first  write  (21)  In  terms  of  three  first- 
order  equations  by  Introducing  new  dependent  vari- 
ables u(x,n.t),  v(x,n,t).  that  Is, 


(bv)-  . J fv  * m3  . x(u  V ft  ♦ (24c) 

We  next  consider  the  net  cube  showir  In  fig.  1 and 
write  difference  approximations  to  (24)  as  discussed 
in  several  references  (see,  for  example.  Ref.  (9)). 


Fig.  1 Net  cube  for  the  standard  Box  method 

Equations  (24a)  and  (24b)  are  approximated  using 
centered  difference  quotients  and  averaged  about  the 
midpoint  (x^  ,tn,nj.v,).  The  difference  quotients 
which  are  to  approximate  (24c)  are  written  about  the 
midpoint  (x^.ij,  t^.i^,  nj-',)  of  the  cube  whose  mesh 
widths  are  r<,kn  and  )ij.  The  resulting  algebraic 
system,  which  Is  nonlinear.  Is  linearized  by  Newton's 
method,  and  the  linear  system  Is  solved  by  the  block 
elimination  method  discussed  In  Ref.  (9). 

Modified  Box  Method 

When  there  Is  flow  reversal  across  the  boundary 
layer  at  some  x and  t,  we  modify  the  standard 
Box  method  used  for  (24c)  but  retain  that  for  (24a) 
and  (24b)  and  center  them  at  (x< ,tn,nj.),).  To  write 
the  difference  aporoxlmatlons  for  the  Box  centered  at 
(xi-ij,  tn-i,,  '’i-','*  ***  examine  previously  computed 

values  of  uj’ij.  If  * 0»  then  we  use  the 

standard  Box  method:  If  uT:!!  < 0,  then  we  write 
(24c)  for  the  Box  centered  at  P (see  Fig.  2)  using 
quantities  centered  at  P,  Q,  and  R,  where 

**  * ^*1’  ‘n-S* * ^*i-S,*  ^n*1j-S^’ 

R * ^*1+i,*  *n-r 

Equation  (24c)  can  then  be  written  as 


Fig.  2 Finite  difference  molecule  for  the  zig-zag 
differencing 

(bv)'(P)  ♦ (fv)(P)  ■ x(P)  9u(Q)  (Q)  ♦ ♦u(R)  |^(R) 

-ov(Q)  (Q)  - ♦v(R)  (R)  ♦ 1^  (P)|  (26) 
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The  resulting  algebraic  system  Is  again  nonlinear  and 
its  solution  Is  obtained  by  using  the  procedure  fol- 
lowed in  the  standard  Box  method. 

RESULTS  AND  DISCUSSION 

To  study  the  calculation  of  unsteady  turbulent 
flows  with  flow  reversal,  we  chose  an  external  veloc- 
ity distribution  of  the  form 


lie  • %(!  - <>(x  - x‘)(t  - f’))  0 < X < 1.  t>0 

(28) 

where  a Is  a positive  constant.  The  same  velocity 
distribution  was  recently  used  by  Cebeci  (3)  for  lam- 
inar flows  in  order  to  study  the  computation  of 
unsteady  laminar  flows  with  flow  reversal  using  the 
solution  procedure  described  earlier,  and  to  see 
whether  there  Is  a singularity  associated  with  such 
flows.  Figure  3 shows  the  computed  local  skin-friction 
coefficient  and  Fig.  4 shows  the  corresponding  results 
for  displacement  thickness  In  0 < x < 1 for  various 
values  of  t with  a • 20.  It  Is  clear  from  these 
results  that  the  solution  remains  smooth  for  all  cal- 
culated values  of  t,  even  when  the  region  of  rever- 
sed flow  occupies  the  majority  of  the  boundary  layer 
and  there  Is  no  hint  of  sinqularlty.  Instead,  as  In 
the  study  reported  In  Ref.  (S),  we  see  the  familiar 
rapid  Proudman-Johnson  qrowtF  (10)  of  the  boundary 
layer  In  the  reversed  flow  region  which,  unless 
special  measures  are  taken,  terminates  the  calculations 
due  to  the  rapid  thickening  of  the  boundary  layer. 

Figures  5 and  6 show  the  calculations  for  turbu- 
lent flows  with  a • 40  and  a unit  Reynolds 
Ug/v  • 2.2  X 10*/m.  The  results  shown  In  Fig.  5 were 
obtained  by  using  different  expressions  for  A;  those 
shown  by  circles  were  obtained  by  (8),  and  those  shown 
by  solid  lines  by  writing  (8)  as 
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Fig.  4 Variation  of  dimensionless  displacement 
thickness  «*/L/Rl  various 

values  of  t. 


A •26(1);''  (29) 

As  in  seen,  both  expressions  give  nearly  the  same  results . 

The  results  In  Fig.  6,  as  in  laminar  flows, 
indicate  no  signs  of  singularity  for  all  calculateil 
values  of  t.  This  is  In  contrast  to  the  findings  of 
Patel  and  Nash  (11).  Again  we  see  the  familiar  rapiii 
thickening  of  the~boundary  layer  in  the  reversed  flow 
region.  If  It  were  not  for  this,  the  calculations  for 
greater  values  of  t than  those  considered  here, 
would  have  been  computed. 

Table  1 shows  a comparison  between  Cf-values 
obtained  by  CS  and  BF  models.  At  the  present  the 
calculations  with  the  BF  model  have  only  been  done  for 
positive  wall  shear.  For  this  reason,  the  comparisons 
In  Table  1 are  limited  to  those  for  which  Is 
positive. 

As  seen  from  the  values  In  Table  1,  the  predic- 
tions of  the  CS  and  the  BF  models  are  nearly  the  same 
for  steady  flow.  Initially  the  Cf  value  at  x • 0 
calculated  by  the  CS  model  Is  slightly  lower  than  that 
computed  by  the  BF  model;  this  Is  due  to  matching. 
However,  at  some  distance  downstream,  both  methods 
give  nearly  the  same  answers.  For  unsteady  flows. 


Mg.  S Local  sk  in-f  net  ion  variation  with  x for 

various  values  of  t.  Solid  lines  denote  the 
calculations  made  by  (29)  and  circles  by  (B). 


Fig.  6 Variation  of  displacement  thickness  with  x 
for  various  values  of  t. 

with  pressure  gradient  Increasing  as  a function  of 
time  and  distance,  the  differences  increase.  Since 
the  same  numerical  method  is  used  in  both  cases,  the 
difference  may  be  attributed  to  the  turbulence  models. 
Additional  studies  are  required  to  explain  the  pre- 
dictions of  each  model  but  the  differences  may  be  due 
to  the  better  ability  of  the  transport  model  to  cope 
with  strong-pressure  gradient  flows. 
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T*ble  1 Computed  c.-v*lues  (x10>)  with 
CS  and  6F  mo^Is 

t*0  t»0.10  t»  0.225 

CS  BF  CS  BF  CS  BF 


0 0.»1  0.318 

0.25  0.292  0.308 
0.50  0.285  0.292 
0.75  0.279  0.279 
0.95  0.274  0.272 

references 


0.301  0.318 
0.275  0.306 
0.248  0.283 
0.250  0.271 


0.301  0.318 
0.214  0.297 
0.102  0,222 
0.167  0.207 


0.277  0.276  0.319  0.302 
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